Structural approximations to positive maps and entanglement breaking channels 
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Structural approximations to positive, but not completely positive maps are approximate physical 
realizations of these non-physical maps. They find applications in the design of direct entanglement 
detection methods. We show that many of these approximations, in the relevant case of optimal 
positive maps, define an entanglement breaking channel and, consequently, can be implemented via 
a measurement and state-preparation protocol. We also show how our findings can be useful for the 
design of better and simpler direct entanglement detection methods. 
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I. INTRODUCTION 

Entanglement is one of the most important, and 
presumably necessary, ingredients of quantum infor- 
mation processing For this reason there is a con- 
siderable interest both in theory and experiments in 
designing feasible and efficient ways of entanglement 
detection. Indeed, there has been a lot of progress 
in this problem recently. The most frequently used 
and investigated entanglement detection methods in- 
clude: i) tomography of the quantum state with lo- 
cal measurements, useful for low dimensional systems 
provided entanglement criteria for the states in ques- 
tion are known @, 0, 0], but impractical for higher 
dimensional systems; ii) methods based on detecting 
only some elements of the density matrix for a contin- 
uous family of measuring devices settings, such as the 
method of entanglement visibility 0; iv) tests of gen- 
eralized Bell inequalities [fj], although there are states 
that despite being entangled do not violate any Bell 
inequality 0, [|| nor any known Bell inequality [§] ; v) 
entanglement witnesses [l(J EH]; vi) direct entangle- 
ment detection schemes, for pure [l2j or mixed states 
and, in particular , using structural approxima- 
tions to positive maps flaj; vii) "nonlinear" entan- 
glement witnesses [16| and viii) methods employing 
measurements of varianc es Il7l| or even higher order 
correlation functions [J, Il8|, or relying on entropic 
uncertainty relations [19J. The methods v) and vi) 
are the subject of the present paper and we discuss 
them in more detail below. First we recall some basic 
definitions. 

a. Entanglement Witnesses. An observable 
E = is called an entanglement witness if and only 
if, for all separable states cr, the average tr{Eo) > 
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and there exists an entangled state g for which 
tr(Eg) < 0. As shown in Ref. (lfjfl . the Hahn-Banach 
theorem implies that for every entangled state g, there 
exists a witness E that detects it, i.e. tr(Eg) < 0. 
Conversely, the state a is separable if and only if for 
all witnesses it holds tv(Ea) > 0. As has been pointed 
out in Ref. [131, entanglement witnesses can be effi- 
ciently measured with local measurements and, more 
importantly, one can optimize the complexity of this 
measurement with respect to, for instance, the num- 
ber of measuring device settings. Nowadays, entangle- 
ment witnesses are routinely used in experiments to 
detect entanglement in bipartite [HI and multipartite 
[IHl systems. 

b. Positive Maps. A related concept is that of 
a positive map. Let B(Ha) and B(Hb) denote the 
spaces of bounded operators on Hilbert spaces Ha 
and Hb respectively. Then a linear map A: B(Ha) — * 
B(Hb) is called positive if A(g) > for every g > 0. 
However, not every positive map can be regarded as 
physical, describing e.g. a quantum channel or the 
reduced dynamics of an open system: a stronger pos- 
itivity condition is required [23]. Namely, a map A 
is physical whenever it is completely positive, which 
means that the extended map 1 ® A: B(/C (g) Ha) — * 
B{fC ® Hb) is positive for any extension K. 

Again, as shown in Ref. [loj (see also (HI), a 
state g £ B(Ha ®Hb) is entangled if and only if 
there exists a positive, not completely positive map 
A: B(H B ) -> B(Ha) that detects g, i.e. [1 <g> A](g) is 
not positive definite. A paradigm example of a posi- 
tive but not completely positive map is transposition, 
T, whose great significance for separability was first 
realized in Ref. [21|- It turns out to detect all the 
entangled states in B{C 2 <g> C 2 )and S(C 2 <E> C 3 ) [13]. 
However, as it is well known [26| (see also e.g. Ref. 
[H] and references therein) , in higher dimensions there 
are entangled states which possess the positive partial 
transpose (PPT) property. 
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Entanglement witnesses and positive maps [22j are 
related through the Jamiolkowski isomorphism [28] . 
Let E € B(Ha <8> Hb)- From this moment on 
we assume that the considered Hilbert spaces are 
finite dimensional, dimHA.B — <Ia,b < oo (for 
an example of infinite dimensional generalization of 
Jamiolkowski isomorphism see Ref. [HI])- We then 
define Ae : B{Ha) — ► B(Hb) as follows: 



A E {g) = d A tr A [E(g T ® 1)] . 



(1) 



Conversely, introducing a maximally entangled vector 
in Ha <8> Ha, 



i=X», p + = !*+><*+! (2) 



1$+) = 



we define for each map A: B{Ha) —> B(H E ) an oper- 
ator 



E A = 1 ® A(P + 



(3) 



acting on B(Ha ® Hb)- Then A is positive if and 
only if E\ is an entanglement witness, since E\ E = E 
[28]. Moreover, A is completely positive if and only if 
E\ > 0, i.e. E\ is a (possibly unnormalized) state. 

c. Structural Physical Approximation. 
Positive maps are stronger detectors of entanglement 
than the corresponding witnesses, despite the de- 
scribed Jamiolkowski isomorphism. They detect the 
same states as the corresponding witnesses plus those 
obtained through local invertible transformation on 
one side. Unfortunately generic positive maps are 
not physical and their action cannot be directly 
implemented. It is therefore challenging to try to 
find a physical way to approximate the action of 
a positive map. This is the goal of the structural 
physical approximation (T3. fl5|. The idea is to mix a 
positive map A with some simple completely positive 
map (CPM), making the mixture A completely pos- 
itive. The resulting map can then be realized in the 
laboratory and its action characterizes entanglement 
of the states detected by A. In the particular example 
studied in Refs. [3, [H[ this idea has been applied to 
the map A = 1 ® T. Although experimentally viable, 
this method is not easy to implement since, at least 
in its original version, it requires highly nonlocal 
measurements: subsequent applications of A followed 
by optimal spectrum estimation. A more detailed 
discussion on this entanglement detection scheme is 
given m Section EO 

In the case of finite dimensional Hilbert spaces, we 
can, without lost of generality, restrict our attention to 
contractive structural approximations, i.e. tr(A(p)) < 
1, for tr(g) = 1. If the initial map is not contractive, 
we can always define A'(g) — A(g)/tr(A(g*)), where 
tr(A(g*)) — max e:tr ( e ) =1 tr(A(g))) and the maximum 
is attained for some g* due to the compactness of the 
set of all states. Contractive CPM's can be realized 
probabilistically as a partial result of a generalized 
measurement (see Ref. [HI]). 



d. Entanglement Breaking Channels. A dif- 
ferent use of maps in the context of quantum theory 
concerns the description of quantum channels, which 
are completely positive and trace-preserving maps. A 
channel A that transforms any state g into a separable 
state A{g) is called entanglement breaking (EB) [30]. 
Clearly, these channels are useless for entanglement 
distribution. In Ref. [30], the following equivalence 
was obtained: 

1. The channel A is EB. 

2. The corresponding state E\ is separable. 

3. The channel can be represented in the Holevo 
form: 



A(g) = ^t4F k g)g k , 



(4) 



with some positive operators F k > defining a 
generalized measurement [3l|, J\ F k — 1, and 
states g k determined only by A. 

The last property above means that the action of an 
EB channel can be substituted by a measurement and 
state-preparation protocol. Moreover, from the sepa- 
rable decomposition of the state E\ 



(5) 



with \v k ) £ Ha and \w k ) £ Hb, one obtains the fol- 
lowing explicit Holevo representation of A [3(3] : 

A(g) = ^\w k )(w k \tT {d A p k \v k ){v k \) g , (6) 



where the overbar denotes the complex conjuga- 
tion. The positive operators {dAPk\v k )(v k \} define a 
properly normalized measurement due to the trace- 
preserving property of A. 

Notice that these results can easily be extended to 
the case of contracting maps. Then, a Holevo de- 
composition is still possible for EB maps with the 
positive operators F k defining a partial measurement, 

In this paper we address the question of imple- 
mentation of structural approximations to positive 
maps through (generalized) measurements. In par- 
ticular, we study structural approximations to maps 
A: B(Ha) — * B(Hb) obtained through minimal ad- 
mixing of white noise: 



A(p)=ptr{p)— + (l-p)A(p). 



(7) 



Minimal means here that we take the smallest noise 
probability < p < 1 for which A becomes completely 
positive. Now the key question is when such A can 
be implemented through generalized measurements, 
i.e. when they correspond to EB maps according to 
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Eq. |(4|). As a consequence, we are led to study the 
separability of witnesses of the form 

E A = 1® A(P+) = -P—l + {l-p)E K (8) 
uaclb 

for minimal p such that 

Ea > 0. (9) 

Recall that this is equivalent to A being completely 
positive. In general, we will consider contractive maps 
A and ask whether they correspond to EB maps, not 
necessarily trace-preserving. Note that if a CPM A 
is contractive and EB, then there exists an EB ex- 
tension to a trace-preserving map, A'(p) = A(g) + 
[tr(p) — tr(A(g))]l/d. In the language of witnesses, 
trace preservation means IvbE^ = 1a- 

The main subject of the paper is the following 
conjecture: 

Conjecture: Structural physical approximations to 
optimal positive maps correspond to entanglement 
breaking maps. Equivalently, structural physical 
approximations to optimal entanglement witnesses E 
are given by (possibly unnormalized) separable states. 

We prove the above conjecture in several special 
cases and discuss a large number of generic examples 
providing evidence for its validity. This is done for 
both decomposable entanglement witnesses, that de- 
tect only entangled states with negative partial trans- 
position (NPT), and also for non-decomposable wit- 
nesses which in addition detect PPT entangled states. 
Once more, we expect the conjecture to be valid 
in general, but in some cases we restrict our exam- 
ples to structural approximations which are trace- 
preserving. Note that such restriction makes the con- 
jecture weaker, since every contractive EB channel has 
a trace preserving EB extension, but not vice-versa. 

The importance of our result is twofold: i) if the 
conjecture is true, structural physical approximations 
to optimal maps admit a particularly simple exper- 
imental realization — they correspond to generalized 
measurements (33 |; ii) the results shed light on the 
geometry of the set of entangled and separable states 
(cf. Ref. [H). 

The paper is organized as follows. In Section 
Ull we recall the notions of decomposable and non- 
decomposable entanglement witnesses and their op- 
timality, based on the Refs. [sHODl]. In Section HTT1 we 
concentrate on decomposable maps. First we study 
dimensions 2® 2 and 2® 3, where the positivity of the 
partial transpose provides the separability criterion. 
Here we show that in general, without the assump- 
tion of optimality, the conjecture is not true, while it 
obviously holds for optimal decomposable witnesses. 
Next, we discuss general decomposable maps in 2 ® 4 
systems, which are nontrivial due to the existence of 
PPT entangled states. Other examples of maps in 
3 ® 3 systems satisfying the conjecture are presented 



in Appendix B. We conclude this section proving the 
conjecture for the transposition and reduction map 
[38] in arbitrary dimension. Section [IV] is devoted to 
non-decomposable positive maps. We start the dis- 
cussion by analyzing the case of Choi's map, one of 
the first examples of a map in this class. Then, we 
study a positive map based on unextendible product 
bases (UPB's) (36j. Finally, we end this section with 
an analysis of the Breuer-Hall map [39L Hjj, which 
can be understood as the non-decomposable version 
of the reduction criterion. Here symmetry methods 
turn out to be indispensable. We introduce and study 
in some detail a new family of states — unitary sym- 
plectic invariant states. The most technical details of 
these states are mainly given in Appendix C, where, 
as a byproduct, we show that this family includes also 
bound entangled states. Finally, we study the physical 
approximation to partial transposition, as this map is 
used in the direct entanglement detection method pro- 
posed in jlji]. In the latter case the analysis is again 
made possible due to symmetry arguments, in partic- 
ular the unitary UUVV symmetry (cf. Refs. (ill. 42]). 
The paper ends with the conclusions in Section [VTl 



II. OPTIMALITY OF POSITIVE MAPS AND 
ENTANGLEMENT WITNESSES 

The notion of optimality of positive maps and en- 
tanglement witnesses has been introduced in Refs. 
[34L I35I ]. We review it here without proofs, which 
can be found in the original papers. There are two 
concepts of optimality: one general, and one strictly 
related to non-decomposable positive maps (or entan- 
glement witnesses) and PPT entangled states. We 
focus below on entanglement witnesses — the trans- 
lation to positive maps is straightforward using the 
Jamiolkowski's isomorphism (cf. Eqs. HJ and ([3])). 



A. General Optimality 

Let us introduce the notion of general optimality 
first. Given an entanglement witness E we define: 

• D E — {g > : tr(Eg) < 0} — the set of opera- 
tors detected by E. 

• Finer witness — given two witnesses E\ and E2 
we say that E2 is finer than E\, if De x C De 2 , 
i.e. if all the operators detected by Ex are also 
detected by E2. 

• Optimal witness — E is optimal if there exists 
no other witness which is finer than E. 

• Pe = {u®v e Ua®Ub ■ (u®v\Eu®v) = 0} — 
the set of product vectors on which E vanishes. 
As we will show, these vectors are closely related 
to the optimality property. 
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Vectors in Pe play an important role regarding entan- 
glement. A full characterization of optimal witnesses 
is provided by the following theorem: 

Theorem 1: A witness E is optimal if and only 
if for all operators P > and numbers e > 0, E' = 
E — eP is not an entanglement witness. 

In this paper we will use the following important 
corollary: 

Corollary 2: If the set Pe spans the whole Hilbert 
space Ha <£> Wb, then E is optimal. 

B. Decomposable witnesses 

There exists a class of entanglement witnesses which 
is very simple to characterize — decomposable entan- 
glement witnesses [HI]. Those are the witnesses which 
can be written in the form: 

E = Q 1 + Q T 2l (10) 

where Q\ y 2 > and T refers to partial transposition 
with respect to the second subsystem: 

Q T = (1®T)Q. (11) 

As it is well known, these witnesses cannot detect PPT 
entangled states. We recall here some simple proper- 
ties of optimal decomposable entanglement witnesses: 
Theorem 2: Let E be a decomposable witness. If 
E is optimal then it can be written as E = Q r , where 
Q > contains no product vector in its range. 
This result can be slightly generalized as follows: 
Theorem 2': Let E be a decomposable witness. If 
E is optimal then it can be written as E — Q T , where 
Q > and there is no operator P in the range of Q 
such that P r > 0. 



C. Non-decomposable witnesses 

Entanglement witnesses which are able to detect 
PPT entangled states cannot be written in the form 
(flOj) [24], and are therefore called non-decomposable. 
The present Subsection is devoted to this kind of wit- 
nesses. The importance of non-decomposable wit- 
nesses for detecting PPT entanglement is reflected by 
the following: 

Theorem 3: An entanglement witness is non- 
decomposable if and only if it detects some PPT en- 
tangled state. 

We now recall some definitions which are parallel to 
those provided previously. Given a non-decomposable 
witness E, we define: 

• d E = {g > : g T > and ti (Eg) < 0} — the 
set of PPT operators detected by E. 

• Finer non-decomposable witness — given two 
non-decomposable witnesses E\ and E% we say 
that Ei is nd-finer than E\ if ds 1 C c?e 2 , i.e. if 



all PPT operators detected by E\ are also de- 
tected by E 2 . 

• Optimal non-decomposable witness — E is opti- 
mal non-decomposable if there exists no other 
non-decomposable witness which is nd-finer 
than E. 

Again, vectors in Pe play an important role regard- 
ing PPT entangled states. The full characterization 
of optimal non-decomposable witnesses is given by an 
analog of Theorem 1: 

Theorem 4: A non-decomposable entanglement 
witness E is optimal if and only if for all decompos- 
able operators D and e > 0, E' = E — eD is not an 
entanglement witness. 

Note that in principle non-decomposable optimality 
requires the witness to be finer with regard to PPT 
entangled states only, so that a non-decomposable op- 
timal witness does not have to be optimal in the sense 
of Section III A I However, this is not the case since we 
have the following: 

Theorem 5: E is an optimal non-decomposable 
entanglement witness if and only if both E and E T 
are optimal witnesses. 

and 

Corollary 6: E is an optimal non-decomposable 
witness if and only if E T is an optimal non- 
decomposable witness. 

In Ref. [34| optimality conditions have been derived 
and investigated for the case of 2 ® A^-dimensional 
Hilbert spaces. These conditions are, however, very 
complex and for the purpose of the present work we 
will use Corollary 2 to check optimality, even though 
it provides only a sufficient condition. 



III. DECOMPOSABLE MAPS 

This section is devoted to the study of the conjec- 
ture for decomposable maps. We start by proving the 
conjecture for low dimensional systems, namely 2 ® 2 
and 2 ® 3. Then, we provide some rather general re- 
sults for 2 ® 4 systems. Moreover, Appendix B con- 
tains several relevant examples of decomposable maps 
in 3 3 systems where the conjecture also holds. Fi- 
nally, we prove the conjecture in arbitrary dimension 
for two of the most important examples of decompos- 
able maps, the transposition and reduction maps. 



A. 2 ® 2 and 2 <g> 3 

We begin with general examples in the lowest non- 
trivial dimensions. Take A to be a positive map from 
B(C 2 ) to B(C 2 ) or to B(C 3 ). Recall that ev- 
ery such map is decomposable, i.e. is of the form 
A = Af p +To A2 and that its corresponding entan- 
glement witness can be written as l|10p . We will first 
show that not every structural approximation to A is 
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entanglement breaking. In other words, the optimal- 
ly of the positive map is essential for the conjecture. 
For definiteness' sake we analyze the 2<g)2-dimensional 
case, but the argument also holds in 2 ® 3 systems. 

Let us consider the entanglement witness E\ corre- 
sponding to A and, Qi and Q2 in (fTUj) to be rank-one 
operators of the form: 



" ft 








ft 




- 
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Q-2 



with real positive a and b. Then the witness 



(12) 



Qi 



Ql 



a 

b b+a 

b + a b 

a 



(13) 



is not positive and therefore A is not completely pos- 
itive. From the general form $8$ we obtain that: 



Ek = 



(i-p)(Qi + Q£) 

a' + c 

b' + c b' + a' 

b' + a' b' + c 

a' + c 

a'= (l-p)a, b'=(l-p)b, c=| 

This operator is positive for 

4ft 



P> 



4ft + 1' 



(14) 



(15) 



which is the condition for the structural approxima- 
tion. 

In order to study separability, it is enough to check 
the PPT condition, as it is both a necessary and suffi- 
cient condition in the lowest dimensions [25]. Apply- 
ing partial transposition we obtain that the state (fM)) 
is not PPT, and hence entangled, for 



P< 



4b 



46 + 1 



(16) 



Taking b > a the condition 1|T5]) and lfT6|) can be simul- 
taneously satisfied, thus giving a structural approxi- 
mation which is not entanglement breaking. 

Above we have considered a general positive map 
fromS(C 2 ) to£(C 2 ) orto£(C 3 ), i.e. A = Af p + To 
Afr p . Let us now consider an optimal one: 



A = ToA CP , E A = Q l 



(17) 



It immediately follows that any structural approxima- 
tion to such A is entanglement breaking since 



F T - P -\ 

^ 4 



(l-p)Q>0, 



(18) 



so that Ea is separable. Thus, in the lowest dimen- 
sions any structural approximation to an optimal map 
is entanglement breaking. 



More generally, for arbitrary dimension we immedi- 
ately obtain from Eq. (fl8|) that any structural ap- 
proximation to an optimal decomposable map (fT7j) 
(cf. Section III A[) gives rise to a PPT state. How- 
ever, in principle not necessarily that state is separa- 
ble, i.e. not necessarily A is entanglement breaking, as 
PPT condition is no longer sufficient for separability 
in higher dimensions. 



B. 



2®4 



We now study optimal decomposable maps in 2 ® 4 
dimensional systems. The main characterization of 
such witnesses/maps is given by Theorems 2 and 2' 
from Section III Bl and we will use it extensively in 
what follows. In some cases we will present general 
results, while in others we consider what seem generic 
examples, giving evidence supporting our conjecture. 
Other examples of witnesses in 3 <g> 3 systems fulfilling 
the conjecture are presented in Appendix B. 

Let us then consider systems with dimension 2 ® 4. 
There are only three possibilities in this case, depend- 
ing on the rank r(Q) of the operator Q (cf. Theorem 
2, Section III Bp : r(Q) = 1,2, or 3. Higher ranks are 
not possible as then Q would have a product vector 
in its range and hence the witness Q r would not be 
optimal [34j]. 

When r(Q) — 1, then Q is effectively supported in 
a 2 (g) 2 subspace and the results of Section IIII Al im- 
ply that its structural approximation is entanglement 
breaking. 

When r{Q) = 2 there are two further possibilities: 
Q is supported either in a 2 (g) 3 subspace or in the full 
2 ® 4 space. The first case is again covered by Section 
IIII Al In the latter case, Q can be written as a sum of 
projectors: 



where 



Q — P^p + P x , 



^=|0)|/ 1 ) + |l)|/ 2 ), 

x = |o)|/ 3 ) + |i)|/ 4 ). 



(19) 



(20) 
(21) 



Here |0), |1) is the standard basis in C 2 and fi,...,fi 
are vectors in C 4 . In the most general case of con- 
tractive maps, fi and f% are orthogonal to f% and f^, 
and this consists of the only condition required for 
the proof. Note that if the vectors fi are mutually 
orthonormal, the map is trace-preserving. Projectors 
in Eq. lfT9|) define a decomposition of C 4 into a direct 
sum C 2 © C 2 and, hence, Q has a block-diagonal form 
resulting from a split 2 ® (2 8 2) = (2 ® 2) (2 ® 2). 
Applying the results of Section IIII Al to each of the 
2 (g> 2-blocks we obtain the result. 

We are left with the most interesting case: r(Q) = 
3. Take P a projector on the kernel of Q. The state 
P has rank 5, is PPT and possibly entangled. In this 
case we cannot prove the conjecture in general, and we 
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consider an example where the range of P is spanned 
by the product vectors (1, a) <8> (1, a, a 2 , a 3 ), for all 
complex numbers a. This can be always achieved ap- 
plying a local invertible transformation on C 4 side (cf. 
[37||). Since, by construction, Q is supported on P's 
kernel, it is supported on a span of the vectors 



|10)-|01), 
|02) 

|03)-|12) 



(22) 
(23) 
(24) 



where \ij) denotes the standard product basis of 
C 2 (g) C 4 . As evidence to support our conjecture, 
one can see that the structural approximation to Q T , 
where Q is given by the sum of projectors on the above 
vectors, is indeed entanglement breaking. The details 
of the separability proof are given in Appendix A. 
Note, that the EB map corresponding to Q r is not 
trace-preserving, but as we mentioned in the intro- 
duction it has an EB trace-preserving extension. 



C. Transposition 

We conclude the study of decomposable maps by 
proving the conjecture in arbitrary dimension for two 
of the best known positive maps, the transposition 
and reduction maps. Let us first consider structural 
approximations to transposition T: B(H) — > B(H), 
Ti. = C d , which is an optimal decomposable map, for 
arbitrary dimension. The corresponding witness E T , 
obtained from Eq. <(Sj) , turns out to be a Werner state 
on H <g> H: 



p_. 

d 2 ' 



i-p. 



(25) 



Here F is the flip operator, such that Fip <g) </> = </> <g) 
ip. One easily sees that Et is positive, and hence T 
completely positive, when 



P> 



d+1 



(26) 



which is the condition for the structural approxima- 
tion for T. To check the separability of Et, we use the 
fact that the PPT criterion is necessary and sufficient 
for Werner states. Since, for all p, we have 



P 



E T = ^l + (l-p)P + >0 



(27) 



Et becomes separable at the point it becomes a state. 
This implies that the structural approximation to 
transposition is always entanglement breaking. 

Employing the [/[/-invariance of Werner states, we 
find an explicit expression for T in the Holevo form 
(0]) . Recall that each Werner state can be represented 
using the UU depolarizing map T>uu as [3]: 

Qw = T> vu {q) = [ dU(U ® U)g(rf ® U r ), (28) 



where dll corresponds to the Haar measure over the 
unitary group. Since Werner states are spanned by 
the operators {1, F} [7J, normalized Werner states are 
completely defined by the parameter (F) = tr(gwF). 

For the critical witness, i.e. Et with minimal p, we 
have (F) — 1. One can easily check that the state 
g = 1 00) (00 1 has the same expectation value, hence 



E T = dU\vu){vu\®\wv){wu\ 



(29) 



with \vu) — \u)u) = U\0). Notice that this expression 
is a continuous version of (0, where the discrete set 
of states {|«fc), |wfe)} is replaced by a continuous set 
{|v[/), |u>u)} and the probability pk is replaced by the 
probability distribution dll. According to Eq. ([6]), T 
can be written as 

T(g) = JdU\wu)(wu\tr[(d\vu)(vu\)Q], (30) 

where we used the invariance of the integral under 
conjugation. This approximation has a clear intuitive 
explanation. Given an unknown state, first one tries 
to estimate it in an optimal way using the covariance 
measurement defined by the infinite set of operators 
{My = d\vu)(vu\}, distributed according to the Haar 
measure. If the measurement outcome corresponding 
to \vjj) is obtained, the state \vu){vu\ T = \wu)(wu\ 
is prepared. Finally, it is important to mention that 
the map defining the depolarization process T>uu can 
also be implemented by the finite set of unitary op- 
erators {pk, Uk} of which in our case leads to a 
measurement with a finite number of outcomes. 



D. Reduction Criterion 

Finally, we consider the (normalized) reduction map 
A/? defined as follows: 



Afl(/>) 



1 



d- 1 



tr(p)l - p 



(31) 



which is also an optimal decomposable map [38]. The 
condition for the structural approximation A R to be 
completely positive reads: 



l®A fl (P+) = E t 



d — V 1 — V , . 

l-^-VP + > 0, (32) 



d 2 (d-l) d-1 
which is immediately equivalent to: 

d 



P 



> 



1 



(33) 



In order to study the separability of En, note that En 
is an isotropic state, i.e. Er is [/[/-invariant. For such 
states the PPT criterion is again both a necessary and 
sufficient condition for separability [H, [4l|. Denote 
by n± the projectors onto the symmetric Sym(7i ® 
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H) and skew-symmetric HAH subspaces respectively. 
With the help of the identities P+ = (l/d)F T and 



II + — n_ one obtains that: 



2d-p(d+ 1) 
d 2 (d- 1) 



(34) 



which is positive for all p. Hence, when En becomes 
positive, it also becomes separable which implies that 
the structural approximation to A# is always entan- 
glement breaking. 

Again we use the invariance properties of E R to 
write Ar in the Holevo form (j4]) . These states belong 
to a space generated by {1, P+} and therefore can be 
completely described through the parameter (P + ) = 
tr(gP+), For the critical witness, the expected value 
is (P+) = and a possible separable decomposition of 
the state reads: 



E R = / dU{U® U)\(j>){(j>\{U ® U) 



with |0) = 1 01). According to Eq. ©, 



(35) 



A fl (e) = J dU\w u ){w u \tT{{d\v u ){v u \) e ) (36) 
where \v v ) = U\0) and \wu) =U\l). 

IV. NON-DECOMPOSABLE MAPS 



By checking the positivity of this state we find that 
the map Ac is completely positive for p > 3/5. 

The entanglement witness Eq is separable since, for 
critical p, it can be represented by the following convex 
combination of (unnormalized) product states: 



En = 



= j- 5 (a 01 



Ci2 + &02 + <?d) 



(39) 



Here a d = |02)(02| + |10)(10| + |21)(21| is obviously 
separable and the matrices Oij are defined on the sub- 
space ij, i.e. spanned by \jj)}, and 
read: 



oy = 1 - \ii){jj\ - \jj)(n\ 



(40) 



We can easily check that these density operators 
are PPT and hence separable. Choi's map is not 
proven to be optimal, and there are even reasons 
to believe that it is not. Namely, if one looks 
at product vectors at which the mean of E c van- 
ishes, they have the form (1, exp(i^i), exp(i(f>2)) 
(1, exp(— i<f>i), exp(— ifo)), and are orthogonal to the 
vector (1, 0, 0, 0, 1, 0, 0, 0, 1), so that they do not span 
the whole Hilbert space and do not fulfill the assump- 
tions of the Corollary 2 from Section 2. Still, as we 
have shown, the structural physical approximation for 
the Choi's map is entanglement breaking. Thus, if this 
map is (is not) optimal, this supports (does not con- 
tradict) our conjecture. 



In this section, we move to non-decomposable maps. 
We first consider the Choi map, which is one of the 
first examples of a non-decomposable positive map. 
After this, we study those maps coming from unex- 
tendible product bases. Finally, we analyze a recently 
introduced positive map, the Breuer-Hall map. In all 
the cases, we are able to prove the conjecture in arbi- 
trary dimension. 

A. Choi's map 

We now move to a non-decomposable map proposed 
by Choi 0]. The normalized map A c : B(C 3 ) -> 
B(C 3 ) can be written as: 

(37) 

where |i) is a fixed basis of C 3 and the summation is 
modulo 3. According to Eq. ([8]), the witness E c asso- 
ciated with the structural approximation Ac reads: 

Ec - p\ + (38) 

+ (e^^i+i*' 4 - 1 )^*- 1 !]- 3 ^) ■ 



B. UPB Map 

Let us now focus on unextendible product basis 
[36] . Recall that an unextendible product basis in 
an arbitrary space Ha ®Hb — < C dA <8> C ds consists 
of a set of n < d^dg orthogonal product states, 
{vi = Xi ® Vi}i = i, such that there is no product 
state orthogonal to them. It is then impossible to 
extend this set into a full product basis. Given an 
unextendible product basis, one can associate a PPT 
bound entangled state 

1 ™ 
eM = ^^(l^-El^l> (41) 

The state is trivially PPT and entangled as there is 
no product state orthogonal to the UPB. 

A (normalized) witness detecting such states can be 
taken in the form [36] 

%=^|i^^ £ 4 (42) 

where e > 0. A map A M : B{H A ) B(H B ) 
corresponding to E^ can be obtained through the 
Jamiolkowski's isomorphism (cf. Eq. {TJ) and is a 
non-decomposable map since the state Q{ v y is PPT. 
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Let us consider the structural approximation A-{ v \ to 



A{ v j. The witness associated with A{„} reads: 
P 



For convenience, we choose U = J. Now, let S € 
Sp(2n, C)P\U (2n) be a unitary symplectic matrix, i.e. 
a complex matrix satisfying: 



E {v} - 



cLa(Ib 
1 



l + (l-p)E {v} 
I np — edj^dg 



(43) 



n — ed,Adi 



d A d L 



l + (l-p)J2\vi){vi\ .and 



Since by definition all vectors Vi are product, = 

1a <8> A{„}(P+) is separable, once it becomes positive. 
Therefore, structural approximations to positive maps 
(|42|) arising from UPB's are entanglement breaking. 

Since is already in a product state form, the 
Holevo representation comes directly from Eq. J6]) for 
each particular unextendible product basis As 
mentioned, this gives the explicit construction of the 
measurement and state preparation protocol approxi- 
mating the map. 



C. Breuer-Hall Map 



In what remains, we study the Breuer-Hall map, re- 
cently introduced in (39l . |4Q | . This positive map can 
be understood as the generalization of the reduction 
criterion to the non-decomposable case. As we show 
next, this map also satisfies the conjecture for any 
dimension. When proving these results, we are natu- 
rally led to the analysis of a new two-parameter family 
of invariant states, that we call unitary symplectic in- 
variant states. 

For even dimension d — 2n > 4, which from this 
moment on we assume, the reduction map i|3ip can 
be yet improved, leading to the (normalized) Breuer- 
Hall map Sill: 

^bh(p) = t^t [tr(p)l - p - Up T U^ 



d-2 

Here U is any skew-symmetric unitary operator, i.e. 
U'U = 1 and U T = —U. The resulting map is 
no longer decomposable and is known to be optimal 
[39]. From the general formula J8|), the entanglement 
witness associated with the structural approximation 
A BH is given by: 



E 



1 



BH 



d-2 
1-p 



d-2p 
~d?~ 



l-(l-p)P+ 



(1®U)F(1®U^) 



Further analysis of Ebh will be again based on sym- 
metry considerations. First of all, we note that since U 
is non-degenerate (|detf7| = 1) and skew-symmetric, 
there exists a basis, known as Darboux basis, in which 
U takes the canonical form: 



1 

-1 



sfls 



SJS 1 = J. 



(47) 



(48) 



Then: 



S<E>SE BH S*®S T 

= al+/3S® SP+S^ ® S T + 7 (5 <8> SJ)F(S^ J f 5 T ) 
= al + (3 P+ + 7 (1 (g) J) (S ® S)F(S^ <g> 5 f ) (l ® J f ) 
= al+ ( 3P++7(l® j)jP(l® J f ) =Ebh, (49) 

where we introduced constants a, /3, 7 for simplicity 
(cf. Eq. ((45])) . In the second step above we used the 
property: 



SJ = JS (50) 

(it follows from Eqs. (|4"7jl. {48]) and the fact that J = 
J), together with the £/[/-invariance of P+ Then 
in the last step we used the C/CZ-invariance of F Q. 

Thus, we have just proven that the witness Ebh, 
associated with the Breuer-Hall map, is invariant un- 
der transformations of the form S®S, with S unitary 
symplectic. Equivalently, its partial transpose Eg H is 
invariant under transformations of the form S®S. We 
will generally call such operators unitary symplectic 
invariant, or more specifically S S- and S'S'-invariant 
respectively. 

To our knowledge, these operators have not been 
studied systematically as an independent family. 
(44) They form a subfamily of 5C/(2)-invariant states of 
Ref. [46J (see also Ref. [3§| where a subfamily of 
SS- invariant states was introduced and Appendix B), 
but since the number of parameters of the latter fam- 
ily increases with the dimensionality, it is manageable 
only for low dimensions. Below we describe unitary 
symplectic invariant states in any even dimension (see 
e.g. Ref. for a general theory of states invari- 

ant under the action of a group G). The results are 
then applied to the investigation of the entanglement 
breaking properties of the structural approximations 
to the Breuer-Hall map. 

(45) 



1. Unitary Symplectic Invariant States 

In the next lines, we characterize the family of Uni- 
tary Symplectic Invariant states. For the sake of clar- 
ity, here we state the main results, the corresponding 
proofs are then presented in Appendix C. 
(46) First of all, one should identify the space of Hermi- 

tian S'S'-invariant operators. As shown in Appendix 
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C, the spaces of SS- invariant and SS-invariant oper- 
ators are j48|: 

SS-im = Spanjl^P; 7 }, (51) 
SS-inv = Span{l,P+,F J }, (52) 

where A J = (1 <g> J)A(1 <g> J f ). 

For a later convenience we introduce two equiva- 
lent sets of generators given by the following minimal 
projectors: 



and 





n = 




(53) 




Hi = 


I(1_F)-P;', 


(54) 




n 2 = 


ld + F), 


(55) 


n 


= n^ 


= P+, 


(56) 


rii 


= u{ 


= 1(1-F J )-P + , 


(57) 


ri 2 


= n 2 J 


= \(1 + F J ). 


(58) 



Relations (|112M13j) imply that both sets define a pro- 
jective resolution of the identity: 

n a u = 6 a0 n and 5^n a = i, (59) 



and analogously for fl a . Moreover [II a ,ri/3] = [491 ]. 

Projectors II Q and li Q form extreme points of the 
convex set of positive unitary symplectic invariant op- 
erators. This allows us to easily describe the con- 
vex sets S and £ of SS- and SS- invariant states 
respectively. The normalization implies that each 
family of states is uniquely determined by two pa- 
rameters: tv(gF), tr(gPf) for S'S'-invariant states and 
tr(pF J ) ] tr(pP+) for SS-invariant ones (compare with 
Refs. [4l|, where orthogonal invariant states 

were characterized). The extreme points of £ and 
£ are given by the normalized projectors n a /tiTI a 
and IT Q /trII a respectively. We stress that both sets 
live in two different subspaces of the big space of all 
Hermitian operators: £ C Span K {l, F, P^} and £ C 
Span R {l. P + , F J }. Partial transposition T brings one 
set into the plane of the other and allows one to study 
PPT and separability. 

Figure Q] shows the plot of £ together with £ r — 
the set of partial transposes of SS-invariant states. 
For definiteness' sake we have chosen to study partial 
transposes of SS-invariant states, but as we will see 
the situation is fully symmetric. The plane of the plot 
is the space of all Hermitian SS-invariant operators 
with unit trace. The set £ r is given by the convex 
hull of the normalized operators li^/trII Q : 



conv 



AS 



trlli trn 2 trn 3 



C Span M {l,F,P^}. 

(60) 



1/d 




\ 




I 


\ 

\PPT BE 


Werner states 




/ BH 





1/d 1 



<F> 



FIG. 1: The plot of the set £ of S'S'-invariant states 
together with E r — the set of partial transposes of SS- 
invariant states. The thick line with the arrow represents 
the partially transposed witness E S g H {p). The dashed line 
represents Werner states; its prolongation to the vertex 
(d, — 1/d) = IIq gives NPT isotropic states. The family 
p(A) from Ref. [39J is given by the edge, connecting ver- 
tices (-1,1) = n and (1,0) ee n 2 . The plot of S, E r , 
and Ebh(p) is identical, with the axes labels changed to 
{F J ) and (P+) respectively. 



The mentioned symmetry between the families mani- 
fests itself in the fact that by changing the axes la- 
bels (F) ->■ (F J ) and (Pf) -> (P+) one obtains 
the plot of £ and £ r — it is given by the identical 
figure in the corresponding plane. This stems from 
the following observations: tr(Il a F J ) — tr(n Q F), 
tr(n Q P + ) = tv(U a P(), trft Q = trll a = trf[^, and 
tx(UlF J ) = tr^F), tr(l£P+) = tr^Pf ). 

The intersection £ r n£ describes those SS-invariant 
states with positive partial transpose. As shown in 
Appendix C, not all of them are separable, i.e. there 
are PPT entangled states in the family. The extreme 
points of the intersection are given by: 



x Q = (0,0), x 1 
x 2 = (1,0), x 3 



1 

dj ' 
d I 



d + 2 1 d + 2 



(61) 



To prove separability of a given point it is enough to 
show that there exists a normalized product vector 
| it) <8> \v) with the identical expectation values of F 
and Pjj^, for the latter values characterize the state 
uniquely. Using this fact, one can see that the ex- 
treme points of the separability region are x , x\ and 
X2- The remaining part of the PPT region contains 
entangled states. 
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2. Entanglement breaking property of Abh 



rable state \ip) = \4>) ® IV")) where 



We can now return to the study of the witness Ebh 
associated with the Breuer-Hall map (cf. Eq. (|45| 
with U — J). As we have shown in Section IIV C| 
Ebh is a S'S-invariant hermitean operator. Before 
analyzing when it becomes positive, note that: 



= _(|0) + |1) + |2) + |3)) 



E BH = 



1 



d-2 



d-2p 1 



j)E BH {l®J f ). 



(62) 



Thus, E B h > if and only if Eg H > 0, i.e. the 
structural-approximated witness is a PPT state. 

From Eqs.([45]l and (|62jl we obtain that when < 
p < 1: 

-1 < tr(E BH F J ) = tv{E T BH F) < i, (63) 
-~ < tr(E BH P + ) = tr(E r BH Pl) < i (64) 

The corresponding interval p i— > Eb H (p) is depicted 
in Fig. [1] by the thick line with the arrow. We have 
plotted -Z% ff (p) rather than Ebh{p)- One sees that 
the line enters the positive region S at the point 
.t = (0,0), that is when both averages ifM)) and ([64)) 
vanish. Equating any of the expectation values to zero 
gives the condition for the structural physical approx- 
imation: 



P 



> 



d+1 



(65) 



Notice that it is the same bound as in Eq. ([33]) for 
the reduction map. Observing Fig. Q] it is clear that 
any structural approximation to Breuer-Hall map is 
entanglement breaking since the positivity region of 
Eg H is inside the separability region of SS-invariant 
states. 

As a byproduct, we also obtain the minimum eigen- 
value Xmin of the witness Ebh , corresponding to the 
original positive map f44|) . From Eq. ((Sj) it follows 
that at the critical probability p = d/(d + 1) one 
must have p/d 2 + (1 — p)X m i n = 0. This leads to 
Xmm = — 1/rf, which corresponds to the eigenvector 
1$+). Note that this eigenvector shares the symmetry 
of E BH : S ® S\$+) = |$+). 

Again, we are able to provide a representation of 
the structural approximation to Breuer-Hall map us- 
ing the S'5-invariance of the corresponding witness: 



E 



BH 



t\S(S®S)\ip){<p\{S®S)\ (66) 



These states are parameterized by (P+) and (F J ) and, 
for the critical witness Ebh we have (P+) = (F J ) = 0. 
The same expected values are obtained by the sepa- 



1 

71 



(10) - |2» 



Then, the Holevo form of Abh is: 

^bh{q) = J dS\ws)(ws\tr(d\vs)(vs\g) 
with \vs) = S\c/)) and \ws) = S^). 



(67) 
(68) 

(69) 



V. ENTANGLEMENT DETECTION VIA 
STRUCTURAL APPROXIMATIONS 

Before concluding, we would like to discuss the ap- 
plication of these ideas to the design of entanglement 
detection methods. Indeed, one of the main motiva- 
tions for the introduction of structural approximations 
[HI was to obtain approximate physical realizations of 
positive maps, which can then be used for experimen- 
tal entanglement detection. 

The original scheme proposed in works as fol- 
lows, see also Fig. [2j Given TV copies of an un- 
known bipartite state, qab, the goal is to determine, 
without resorting to full tomography, whether the 
state is PPT. The idea is to apply the structural 
approximation to partial transposition to this initial 
state and estimate the spectrum (or more precisely, 
the minimal eigenvalue) of the resulting state using 
the optimal measurement for spectrum estimation de- 
scribed in [H2] ■ Note that the structural approxima- 
tion 1 ® T "simply" adds white noise to the ideal oper- 
ator g r . Thus, it is immediate to relate the spectrum 

of (1 ® T){qab) to the positivity of the partial trans- 
position of the initial state. 

Inspired by the previous findings, we study in this 
section whether the structural approximation to par- 
tial transposition defines an entanglement breaking 
channel. This map is of course not even positive 
(so it does not entirely fit with our main considered 
scenario), but obviously by adding sufficient amount 
of noise it can be made not only positive but also 
completely positive. As we show next, the struc- 
tural approximation to partial transposition does in- 
deed define an entanglement breaking channel when- 
ever dA > ds, which includes the most relevant case 
of equal dimension dA = ds- 

This implies that the entanglement detection 
scheme of Fig. [2ja can just be replaced by a sequence 
of single-copy measurements, see Fig. [2jb, being the 
measurement the one associated to the Holevo form 
of the entanglement breaking channel. This alterna- 
tive scheme is much simpler from an implementation 
point of view since it does not require any collective 
measurement, though the measurements are not pro- 
jective. Moreover, it can never be worse than the pre- 
vious method, and most likely is better (see also [32J). 
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(a) 



(b) 



p - SPA(\®T) - 



p - SPA(1®T) - 



p - SPA(\®T) - 



p- M — 



m. 



p— M — 



M -m 



p— M — n h 



becomes: 



(d A d B ) 2 
P 



(d A d B ) 2 
P 

(d A d B ) 2 



d B 

Q AA ' ® n 
i-P 



P 



iAA 



n 



•>AA' 



BB ' 



n 



n 



BB ' 



(72) 



Since only the last term can be negative, one obtains 
the following condition for structural approximation: 



P> 



d\d B 



d\d B 



1 



(73) 



FIG. 2: The original scheme for direct entanglement de- 
tection proposed in fl5| is shown in (a). Given N copies 
of an unknown state g, it consists of, first, the structural 
approximation of partial transposition acting on the ini- 
tial state, followed by optimal estimation of the minimal 
eigenvalue of the resulting state. In the new scheme, all 
this structure is replaced by single-copy measurements on 
the state. The minimal eigenvalue should then be directly 
estimated from the obtained outcomes. 



Comparison of the above threshold with the one given 
by Eq. ((26)1 with d — d B , shows that in order to make 
1a®Tb completely positive one has to add more noise 
than to make the transposition T alone completely 
positive and hence implementable. In other words, 

1a <8> T B is less noisy than 1^ <8> T B . 

We proceed to study the separability of E x ^t- We 
begin by finding the partial transposition of 
with respect to the subsystem A'B' [45|: 



A. Structural Approximations to 1 g> T 



E 



P 



1 



1<8>T 



(d A d B ) 2 



P_ F AA' 



P 



BB' 



(74) 



Let us then consider the structural approximation 
to transposition extended to some arbitrary auxiliary 
space: 1a ®Tb [HI- Note that, unlike in the previous 
cases, the initial Hilbert space describing the system is 
now explicitly a product H = U A ® U B = C dA <g> C dB . 
Moreover, generally 1 ® A is not the same as 1 ® A, 

although 1(g) A(P+) = 1 <g> A(P+), so this problem 
does not reduce to the previous one. Calculating the 

witness corresponding to 1 A <8> Tb one obtains: 



Applying the same technique as above (cf. Eq. (72 
we find that £L T > if and only if: 



P> 



dAd 2 B 
d A d 2 B + 1 



(75) 



Comparing this to the threshold for positivity (f73|) . 
we see that for d A < dB, i.e. when the extension is by 
a space of smaller dimension, there is a gap between 
positivity and PPT. Hence, in this case, for 



{d A P d B y^AA 



1b) ® (1a 
® Ibb' + 



®2i,0] (P. 

l-p pAA' 
d„ r + 



AB,A'B' 
BB' 



. (70) 



where F BB ' is the flip operator onK B ®%= C dB ® 
C dB and p^ B ' A B 5 P AA are projectors onto maxi- 
mally entangled vectors in the corresponding spaces, 



P 



AB.A'B' 



dA ds 



+ 



d A di 



Yl ^2\i3)AB(kl\®\ij)A>B>(kl\. 



i,k=ij,i=\ 



(71) 

The condition for structural approximation, positivity 
of Ei®t, is most easily derived by using the identity 
F BB ' = nf B ' -U BB ' , where Uf B ' is the projector on 
the symmetric subspace Sym(H B ®7i B >), and intro- 



ducing a projector Q+ A 



1 AA >~P^ A '. Then E^ T 



d 2 A d B 



d A d B 



1 



< P < 



d A d% 
d A d 2 B + 1 



(76) 



the witness (|72|l is not separable and the map 1 A (8 T B 
is not entanglement breaking in this region. Recall 
however that this does not represent any counter- 
example to the conjecture as the initial map is not 
even positive. 

In the case d A > ds, we will use symmetry ar- 
guments to prove the separability of Ex®t- From 
Eq. (|70|) it follows that this state is XJUVV- invariant, 
where U E U(d A ), V € U(d B ) (cf. Refs. [HI, 52] 
where C/C/T^y-invariant states were studied). Since 
both groups U(d A ) and U(dn) act independently it is 
easy to convince oneself [41j that the space of UUVV- 
invariant operators is spanned by {1®1, 1<8)F, P+<g>l, 
P + (g) F}. Following the same approach as in sub- 
section IVI[ we prove the separability of Bi®t in the 
AB : A'B' partition by showing that the state can be 
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written as convex sum of product states, i.e. it is has 
the following representation 



dUdV(U A VBUA>V B <)a(U A V B UA>VB< 



vf 



(77) 



(we omit tensor product signs here for brevity) for 
some a separable in the partition AB : A'B 1 . Given 
that states with this invariance are completely de- 
scribed by parameters (1®F), (P+<g>l) and (P + <g)F), 
a must obey the conditions: tr(crl®F) = tr(Pa®Tl<8> 
F) ,J,t(ctP + 1) = tr(E mT P + ®l) and tv(aP + ®F) = 

tr(-Eig)T-P+ ® F). Such state a = \<p){<p\ can be writ- 
ten as 

\<f) = \4>)AB ® 

= (V5oo|00) + V55r|Ql) + Vau|H» 100X78) 



for 



an 



eUd 2 D + 1 



(l + d A ) 



(79) 



1 



d A dl + 1 



(4+d A -rfs(l+^)) (80) 



1 



i A d B 



-(d%+d A ). 



(81) 



Notice that, as expected, a is only well-defined for 

d A > d B - According to Eq. ([6]), the map 1 ® T(g) 
can be written as 



1 T(g) = J dUdV\ 



w uv)( w uv \^(d A d B \vuv ) (vuv \ Q) ■ 



(82) 

where \v uv ) = U ® V\4>) and \w V v) = U ® V\ip). 
Recall also that the integrals over the unitary group 
defining each depolarization protocol can be replaced 
by the finite sums of, e.g., Ref. [431 ] - 

In the case d A = d B = d, we encounter the struc- 
tural approximation to the transposition map ana- 
lyzed in [15]. As mentioned, by providing the rep- 
resentation (|82| we are able to replace the former 
entanglement detection scheme [l5[ by a much less 
resource-demanding one. In the original proposal, n 
copies of T(g) are prepared, followed by optimal es- 
timation of its minimal eigenvalue by means of a col- 
lective projective measurement on the n-copy state. 
Now, one should just perform local measurements in 
the n copies of g with operators defined in ((82]) and 
with that directly estimate the lowest eigenvalue of 
1(g) T. 



VI. CONCLUSIONS 

In this work, we have studied the implementation of 
structural approximations to positive maps via mea- 
surement and state-preparation protocols. Our find- 
ings suggest an intriguing connection between these 



two concepts that we have summarized by conjectur- 
ing that the structural physical approximation of an 
optimal positive map defines an entanglement break- 
ing channel. Of course, the main open question is 
(dis)proving this conjecture. It would also be interest- 
ing to obtain slightly weaker results in the same direc- 
tion, such as proving the conjecture for general opti- 
mal decomposable maps (which seems more plausible 
due to the fact that the conjecture holds for trans- 
position). We have also applied the same ideas to 
the study of physical approximations to partial trans- 
position, which is not a positive map, and discuss the 
implications of our results for entanglement detection. 

We would like to conclude this work by giving a geo- 
metrical representation of our findings (that should be 
interpreted in an approximate way). It is well known 
that the set of quantum states is convex and includes 
the set of separable states, which is also convex, see 
also Fig. [H These two sets are contained in the set 
of Hermitian operators that are positive on product 
states, which is again convex. Entanglement witnesses 
belong to this set. If the conjecture was true, it would 
mean that the set of optimal witnesses would live in 
a region which is "opposite" to the set of separable 
states, in the sense that when mixed with the maxi- 
mally mixed noise, they enter the set of physical states 
via the separability region. 

Finally, let us mention some further open questions. 
It would be interesting to extend our studies and ask 
which classes of positive maps have structural approxi- 
mation that corresponding to partially breaking chan- 
nels (for definition see [53])? Is our conjecture true for 
maps that are not optimal, but atomic jUj], i.e. de- 
tect Schmidt number 2 entanglement (for definition 
see [55j]? What is the relation between optimality, 
extremality (in the sense of convex sets) and atomic 
property? 
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With the help of the positivity condition (J84J) , the 
resulting matrix can be written as: 




FIG. 3: The sets S, Q and W of separable states, quan- 
tum states and operators positive on product states are 
such that S C Q C W. If the conjecture was true, namely 
all structural approximations to optimal positive maps de- 
fined entanglement breaking channels, it would mean that 
optimal positive maps (witnesses) enter the physical re- 
gion, when adding white noise, via the separability region, 
as shown in the figure. 
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where 



— a z \a I 



,(86) 



> 1. 



(87) 



Appendix A: Proof of the conjecture for a 
rank-three optimal witness in 2 (g> 4 systems 

In this appendix, we show that the structural ap- 
proximation to the optimal witness Q T , where Q is the 
projector onto states (|22| . is separable. Following our 
general procedure (cf. Eq. ([8])), the normalized wit- 
ness associated to the structural approximation reads: 
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(83) 



1 + a 
a 



where a = ^f^. 
tive when 



The above operator becomes posi- 



a(l + a) = 1. 



(84) 



To show that at this point the matrix (|83|) becomes 
separable, we first perform a local invertible transfor- 



mation and pass from Q T 
A\ where 



-al to 1® A 



(V 



■al 1 



.4: 



10 

a 

a 

1 



(85) 



Note that since at the critical point l(84|) . a — a > 0, 
it is enough to show that both matrices in the above 
decomposition are separable. The first matrix, which 
we denote by cr, possesses the following continuous 
separable representation: 



2tt 



where 



V#) 



^|V#))(V#)I, 



-l)®(l,e^,e 2 ^,e 3 ^). 



(88) 



(89) 



The second matrix has a (2 ® 2) © (2 ® 2) structure 
with 2 ® 2 blocks being identical and given by 



10 0-1 

10 

k 

-10 1 



(90) 



Since k > 1 the above matrix is PPT and hence sepa- 
rable. Thus, the whole matrix <f86|l is separable, which 
finishes the proof. 



Appendix B: 3 ® 3 systems 

In this appendix, we provide several examples of 
positive maps satisfying the conjecture. Again we 
consider decomposable optimal maps and study case- 
by-case various possible ranks of the Q operator (cf. 
Theorem 2, Section |nB| . 

The case r(Q) = 1, i.e. Q — \tp)(tp\, splits into 
two subcases. When the Schmidt-rank of \ip) is 2, Q 
is supported in a 2 ® 2 subspace and the structural 
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approximation is entanglement breaking by the pre- 
vious results (cf. Section [ill A|) . In the case where 
\ip) is Schmidt-rank 3, we restrict our attention to the 
trace-preserving case, i.e. assume that \ip) is maxi- 
mally entangled. Alternatively, before checking the 
conjecture we apply local transformations and bring 
to the form i.e. we assume that: 



lows: 



1 

71 



|00) 



in) 



122) 



(91) 



Then the corresponding witness E from Eq. Q turns 
out to be a Werner state of dimension d = 3. This 
witness was already studied for arbitrary d in section 
, where we concluded that such structural approxima- 
tion is always entanglement breaking. 

We move to the case r(Q) = 2. Then Q has to be 
supported either in a 2(8 3 subspace or in the full 3(8)3 
space, since in 2(8)2 there is always a product vector in 
every two dimensional subspace and Q would not be 
optimal by Theorem 2 of Section [TlBl The first case, 
when Q is supported in a 2 (g) 3 subspace, is covered 
by Section fill Al In the other case, we do not have a 
general theory, but in a generic case the range of Q is 
spanned by two Schmidt-rank 2 vectors. We can take 
them to be: 



|01) - |10) 
|12)-|21). 



(92) 
(93) 



Obviously, for such a Q it holds Qe (8 e = =>• (e (8 
e\Q r e (8 e) = for any e € C 3 . Since vectors e ® e 
span the whole C 3 (8 C 3 , by Corollary 2 of Section [HA] 
the witness Q T is optimal. 

Again we do not have a general result here, but 
only consider a generic example of Q given by the 
projectors on the above vectors (|92H93p . The normal- 
ized witness corresponding to the structural approxi- 
mation, E\ = ^-^-(q 1 ^ + alj with a = 9 (^ p ) is given, 
modulo the (1 — p)/4 prefactor, by the matrix: 
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It becomes positive at the point a(a 2 
a = V2. 

which gives the critical probability p c — g ^ +1 

To check the separability at the above point Ip5|). 
note that the matrix (|94| can be decomposed as fol- 
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(96) 



The first two matrices are supported in 2 (8 2 sub- 
spaces. Their partial transposes become positive for 
(1 + a) 2 = 1, which is satisfied at the point ([95|) . The 
last matrix is obviously separable. This allows us to 
conclude that the structural approximation (|94| is en- 
tanglement breaking. 

Next, we consider the case r(Q) = 3. Then Q must 
be supported in the whole 3 (8 3 space (otherwise there 
would be a product vector in the range of Q and Q 
would not be optimal by Theorem 2, Section III Bp . 
In lieu of a general theory, we consider a seemingly 
generic example of 



n , 



(97) 



where by n± we denote the projectors onto the sym- 
metric Sym(H C8 H) and skew-symmetric HAH sub- 
spaces respectively. The corresponding normalized 
witness reads: 



l-p 



E A = Q l 



1 

al =2 



(1 +2a)l- 3P- 



(98) 



where a = 



IL 



and we used the identities F — 
- II_ = 1 - 2II_ and F r = dP + . The condition 



for structural approximation, Ea > 0, is equivalent to 

a > 1. (99) 

Note that the structural-approximated witness f98|) is 
an isotropic state of dimension d = 3 and that this was 



15 



already studied for arbitrary d in Sec. lIIIDl There we 
concluded that such witnesses always correspond to 
entanglement-breaking channels. 

We are left with the last case r{Q) = 4. Note that 
generically if we consider P a projector on the ker- 
nel of Q, then r(P) = 5 and the range of P contains 
exactly < 5 product vectors. In general, Q will con- 
tain some product vector in its kernel and therefore is 
not optimal. For this reason, here we consider not a 
generic but a particular Q where optimality is guaran- 
teed by the Corollary 2 of Section III Al We can treat 
C 3 (g> C 3 as a representation space of two spin-1 rep- 
resentations of SU(2). We then consider positive op- 
erators Q supported on a span of the skew-symmetric 
subspace C 3 A C 3 and the singlet [H]: 

* = -^(|02) + |20> - |11)). (100) 



Denoting by J the total spin, Q is supported on 
the sum of J = and J = 1 subspaces, while P 
is supported on the J = 2 subspace. The kernel 
of Q is then spanned by the vectors of the form 
(1, V2a, a 2 )<g>(l, \/2a, a 2 ) for a complex a. By Corol- 
lary 2 of Section III Al Q r is optimal, as vectors 
(1, y/2a, a 2 ) (g) (1, V2a, a 2 ) span whole of the C 3 <g> C 3 . 

As a particular example we consider 



Q = 2II_ + 2P*. (101) 



The structural approximation gives: 
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(102) 



where 



operator: 



9(1 -p) 



(103) 



The matrix (|102|) becomes positive at the point given 
by the conditions (a + 1) 2 — |=0 and I a + 



2 = 0, which is solved by 



a = - 



(104) 



We now prove that at this point the witness (j 102[) 
becomes separable. We consider the partially trans- 
posed witness: 



1 



2(P 7=1 +P J=0 



(105) 



where Pj projects on the subspace of total spin J. 
Using the technique based on the state invariance de- 
scribed in Sec. IIIICI we explicitly construct a sepa- 
rable decomposition for E^. Analogously to the defi- 
nition (|28| . we introduce spin-1 (£> spin-1 depolarizing 



2>{q) = j cW {1) {U) x 

(17) ® {U)\ g [$> (1) (J7) t ® (t/)f|06) 

= ^tr(gPj =2 )Pj=2 + gtr(ePj = i)Pj=a + 

+tr (gP, =0 )P J=0 . (107) 



where & {1 \U) £ 50(3) denotes spin-1 representation 
of U € SU{2). By direct calculation we check that 



0(|O2)(O2|) + 0(|O1)(O1|) 



(108) 



gives, up to a positive constant, the desired operator 
E\. Since separability of E\ is equivalent to separa- 
bility of Ejy, we have thus shown that the structural 
approximation to the map defined by Eq. (|10ip is 
entanglement breaking. 
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Appendix C: Analysis of Unitary Symplectic 
Invariant States 



The scope of this appendix is to provide a charac- 
terization of the properties of SS and SS invariant 
states. The first step is to find the space of Hermitian 
SS-invariant operators. The corresponding space of 
SS-invariant ones is related to the latter by partial 
transposition T. Since unitary symplectic transfor- 
mations S are obviously unitary, all [/[/-invariant op- 
erators are also SS-invariant. As it is well known, the 
former space is spanned by 1 and F Q. As a rule, 
shrinking the group enlarges the space of the invariant 
operators, so one expects more than that. The form 
of the invariance group G = Sp(2n, C)n[/(2n) implies 
that {G - inv} = {Sp(2n, C) - inv} U {U(2n) - inv} 
(in some sense we will not specify here; see Ref. (4ll|). 
Thus, one has to find the Sp(2n, C)-invariant opera- 
tors. 

Let A be Hermitian and such that: 



^ ^ S>ijSkiAji mn S rm S sn A 



(109) 



ikrs • 



for all S from Sp(2n,C) (now S satisfies _Eq. (08]) 
only). Since S and its complex conjugation S are inde- 
pendent for a general S G Sp(2n, C), and the defining 
equation 1|48|1 does not involve complex conjugation, 
the only possibility for Eq. (|109[) to hold is when A 
is rank one, i.e. Aji mn — ipji4> m n- Then Eq. (j 109|) 
becomes: 



T\ - 



4>ik<l>r 



(110) 



But the only quadratic form that S preserves is J, 
which implies that one must have ■f/'ifc = ciJik and 
4>rs = C2J r s for some complex c\.2 ^ 0. We choose 
ci = C2 = — l/\[d, d = 2n, which leads to: 



; , k 



(|10)-|01) + |32)-|23) + ...) 



Vd 
(1® J)|$+ 



(111) 



(cf. Eq. 0). Hence, P{ = (1 ® J)P+(1 ® J f ) is 
the only Sp(2n, C)-invariant operator, up to a mul- 
tiplicative constant (42|. Using this fact we con- 
clude that the space of S'S'-invariant operators is 
spanned by {1, F, P+}. Correspondingly, the space of 
SS-invariant operators is spanned by {1,F,P+} = 
{1, P+, F J } = (1 ® J){1, P{, F J }(1 ® Jt). As a side 
remark, we note that since J is real, J' = J T = — J 
(cf. definition gHJ)) and hence (1 ® J)A(1 ® J f ) = 
-(1(g) J)A(1 ® J) for any A. We will use this fact 
frequently, but keep writing J' , 

As ageneral rule, G-invariant operators form alge- 
bras [4l| . The constituent relations for the algebras of 



FP J 


= -H 


= p(f 




and 




P+F J 


= -P+ 


= F J P 



unitary symplectic invariant operators are as follows: 

(112) 

(113) 

The above relations follow from the identity .F(l ® 
J)|$+) = (J® = -(1 ® J)\$+), equivalent to 

F J |$+> = -|*+). 

Let us now focus on the study of the PPT region, 
resulting from the intersection E r flS. As we men- 
tioned, when studying separability, one should char- 
acterize the expectation value of the generators of the 
group with product vectors. For a vector \u) ® \v) one 
obtains that: 



(F) =\(u\v)\ 



\UqVq + UiVi + 1*2^2 + U 3 V 3 
1 , 



U2nV2n\ , 



1, ,2 
= j \UoVl — UiVo H \- U2n-lV2n — U2nV2n-l\ ■ 

a 1 1 

(114) 

From these equations, one easily sees that the first 
extreme point from l|6ip can be realized by e.g. u = 
(l/2)(-l, 1, 1,1,0,. ..) and v = 1/^2(1,0,0, 1,0, . . .), 
while points X2 , £3 can be obtained from u = 
1/V2(|0> T |1», v = 1/V2(|0) + |1)) respectively. To 
show that only the set convjxo, Xi, x%) is separable we 
will employ the Breuer-Hall map (|44"1) itself. Note that 
the corresponding separable set convjiro, x i, X 2] T C 
SnS r is determined by the points with the same 
coordinates as xq,X\,X2 but in the (P+), (F J )-plane 
(since e.g. tv(g r P+) = 1/d <^>tv(gF) = 1, etc). 

For an arbitrary SS-invariant normalized state g — 
al + f3F + jP( it holds tr B g= (da + F + (l/d)j)l = 
1/d, since trg — d 2 a + d/3 + 7 = 1 and tr^P/ = 
trgPf = 1/d as J is unitary. Analogously, for an 
arbitrary SS-invariant state g = al + (3F J + 7P+, 
trsg = 1/d, since tr^F' 7 = tr^F = 1. Hence, the no- 
detection condition 1 ® h.BH{g) > takes the same 
form for both families: 



-1-Q- (l® J)/(l® J f ) > 0. 



(115) 



We multiply the above inequality by P+ and P+ 
respectively. Noting that Pf,P+ > and [1 ® 
A-bh(q), P+] = = [1 ® A B h{q), P+], we obtain that 
if a state is not detected by the Breuer-Hall map then: 



, l-tr(oF) 
tr(gPi) < ^— L , or 



tv(gP+) < 



1 - tr(gF J ) 



(116) 
(117) 



respectively. Equivalently, states breaking the above 
inequalities, i.e. states lying above the line (Pf) = 
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(1 - (F)) jd, or above the line (P+) = (l - (F J )) jd in 
the case of SS-invariant states, are detected by Abh 
and hence entangled. 

The set of PPT entangled SS-invariant states is de- 
picted in Fig. [TJ Note that when d — > oo, d even, the 
point X3 — > x 2 , cf. Eq. (|6l)l . and the set of PPT bound 
entangled states collapses. Since we expect that away 
from region boundaries in Fig. [T]the properties of SS- 



invariant states are shared by the states in a small ball 
around them, the collapse of the "volume" of the PPT 
states is to be expected according to Ref. jHl|. From 
the previous arguments (cf. remarks after Eq. (|60j) ) 
and Eq. I|117p . the corresponding diagram for S'S'- 
invariant states is identical, modulo the labels of the 
axes. This finishes our analysis of unitary symplectic 
invariant states. 
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